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Understanding the like-sign dimuon charge asymmetry in pp collisions 
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The D0 collaboration has measured the like-sign dimuon charge asymmetry in pp collisions at the 
Fermilab Tevatron collider. The result is significantly different from the standard model expectation 
of CP violation in mixing. In this paper we consider the possible causes of this asymmetry and 
identify one standard model source not considered before. It decreases the discrepancy of the like- 
sign dimuon charge asymmetry with the standard model prediction, although does not eliminate it 
completely. 

PACS numbers: 13. 20. He 



I. INTRODUCTION 

The D0 Collaboration has measured [l|-[|| the like- 
sign dimuon charge asymmetry and the inclusive muon 
charge asymmetry in pp collisions at a center-of-mass en- 
ergy s/s = 1.96 TeV at the Fermilab Tevatron collider. 
After subtracting the known background sources, the in- 
clusive muon charge asymmetry is found to be compati- 
ble with zero while the like-sign dimuon charge asymme- 
try significantly deviates from the standard model (SM) 
prediction. This deviation is usually interpreted as the 
anomalous charge asymmetry of B° and B® semilcptonic 
decays H, @. 

In this paper we consider other possible sources of the 
like-sign dimuon charge asymmetry, taking into account 
the constraint that the inclusive muon asymmetry is con- 
sistent with zero. We identify one significant contribution 
which was not accounted for previously. In addition, we 
discuss all other sources of the dimuon charge asymmetry 
and show that the measurements of the D0 collaboration 
put strong constraints on them. After presenting avail- 
able experimental results related to the dimuon charge 
asymmetry in Section II, we consider in Section III one 
by one the contribution of different processes into this 
asymmetry. Our results are summarised in Section IV. 



II. EXPERIMENTAL RESULTS 

We use the results of Ref. [3| and express them in a 
model independent way as the charge asymmetry of the 
inclusive muon sample as, and the charge asymmetry of 
the like-sign dimuon sample As- We follow the notations 
and definitions of Ref. [4J, where the asymmetries as 
and As are computed using the muons coming from the 
decays of b and c quarks and r leptons, and from decays of 
short-living mesons such as <j), uj, r\, p°, J/ip. The origin 
of the asymmetries as and As may be the semilcptonic 
charge asymmetry of B° and B° s decays, as well as other 
CP-violating processes. 



From the information presented in Q we obtain 

a s = (-0.063 ± 0.079(stat)±0.141(syst))%, (1) 
A s = (-0.383 ±0.092(stat)±0.102(syst))%. (2) 

These values are obtained, respectively, by multiplying 
the results given in Eqs. (34) and (35) by the coefficients 
c b and C b given in Eqs. (30) and (31) of Ref. 0]. An 
important observation which can be derived from these 
results is that as is compatible with zero, while As is 
significantly different from zero. 

The composition of the inclusive muon sample is pre- 
sented, for convenience, in Table fl] taken from Ref. Q. 
The composition of the like-sign dimuon sample can be 
derived from the information given in Table U assuming 
that the two muons come from independent processes. 
This assumption can be applied because the requirement 
that the invariant mass of the two like-sign muons is 
greater than 2.8 GeV, which is used for selecting the 
dimuon pairs suppresses decays with the two muons 
originating from the same b quark. Since the oscillation 
of the D° meson is found to be small H , the contribution 
of process T 6 in Table [Q to the like -sign dimuon sample 
is suppressed and, therefore, is neglected in the following 
discussion. 

For our purposes it is convenient to classify the pro- 
cesses T 3 into three different categories denoted T 3c t(fi), 
T 3s (fj) and T 3 f s . Processes T 3 d(fi) correspond to the de- 
cays of B° or B° to the final state ft containing two 
c quarks and accessible to both B° and B°. At the 
quark level these decays correspond to the process (bd) or 
(bd) —> cedd. Similarly, processes T 3s (fj) correspond to 
the decays of B® or B® to the final state fj containing two 
c quarks and accessible to both B® and B®. These de- 
cays at the quark level are (bs) or (bs) —J- cess. All other 
decays of B hadrons producing two charm quarks are 
flavour specific and included in the group denoted T 3 f s . 
The weights of these processes are respectively w 3 d(fi), 
w 3s(fj) an d W3f s . By definition, 



w 3 



W 3 f s 



(3) 
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TABLE I: Heavy-quark decays contributing to the inclusive 
muon and like-sign dimuon samples taken from Ref. [j. 
The abbreviation "non-osc" stands for "non-oscillating," and 
"osc" for "oscillating." All weights are computed using MC 
simulation, xo = 0.1260 ± 0.0037 is the time-integrated mix- 
ing probability Q- 



asymmetry are 



Process 



Weight 



Ti 

T-2 

T 3 
T 4 
T B 
T 6 



pTX 

->■ fi~X (non-osc) 
-> b — > lT X (osc) 
ti+X 

c — > jj, + X (non-osc) 
-> b — > c — > /i + X (osc) 
ccq with c — > fi + X or c — > /i 

p°,^(1020),J/V,^'-»-/iV 
66cc with c — s> fi + X or c — ► fx X 
cc with c — > fi + X or c — > pT X 



b 
b 

b 

6- 
6 -)• 



Ml = 1. 

Wl a = (1 — Xo)wi 
Wlb = X0™1 

w 2 = 0.096 ± 0.012 

W2a = (1 — X0)W2 
W2b = X0™2 

"X w 3 = 0.064 ± 0.006 
w 4 = 0.021 ± 0.002 
w 5 =0.013 ±0.002 
w 6 = 0.675 ±0.101 



III. CONTRIBUTIONS TO THE LIKE-SIGN 
DIMUON CHARGE ASYMMETRY 

We consider in this section some SM processes pro- 
ducing the like-sign dimuon charge asymmetry As, and 
the current experimental constraints of their contribu- 
tion taken from Q. Particles of physics beyond the SM 
may add new Feynman diagrams with loops in B° -R- B° 
mixing, B® O B® mixing, and in penguin decays. They 
are not discussed here. 

The main source of the like-sign dimuon pairs produced 
in pp collisions is bb events. The background muons from 
K — > llv and 7r — > \iv decays are excluded by definition 
from the asymmetries as and As, and are not discussed 
here. One of the B hadrons from the bb pair decays to 
a "right sign" muon, i.e. to a muon with charge of the 
same sign as the charge of the initial b or b quark at the 
time of the pp collision. To obtain a like-sign dimuon 
event, the other 6-hadron must decay to a "wrong sign" 
muon. 



A. CP violation in mixing of B° and B® mesons 



The complex phase (j>\ 2 °f the mass matrix of the B® 
(q = d, s) system produces the charge asymmetry of 
the "wrong sign" semileptonic B® decays defined as 



-.9 - 



r(B° q (t)->i+x)-T(B° q (t)^rx) 
" sl " r(B°(t) -> i+x) ± r(B»(t) -> i-x) ' 

The asymmetry a q si is related to the phase <j>\ 2 as 



at, = — — - tan ( 



(4) 



(5) 



The contributions of this asymmetry to the inclusive 
muon charge asymmetry and the like-sign dimuon charge 



as(from a q sl ) = c.bC q a% 
4s (from a q sl ) = C b C q a q sl , 



(6) 
(7) 



where the coefficients c&, Cb, Cd and C s are defined in 
Ref. 0: 

c b = 0.061 ±0.007, C b = 0.474 ± 0.032, (8) 

C d = 0.594 ±0.022, C s = 0.406 ± 0.022. (9) 

The SM predictions of the phases 4> q 2 an d asymmetries 



a:, are 



#? 2 (SM) = -0.075 ±0.024, (10) 

<^ 2 (SM) = ±0.0038 ±0.0010, (11) 

«4(SM) = -(4.1 ±0.6) x 10~ 4 , (12) 

<4(SM) = ±(1.9 ±0.3) x 10~ 5 . (13) 



Recently, the experimental measurements of both a^ 
and a*) became available. The measurements of a^, ai 
performed at T(45) [3] and by the DO experiment 
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-0.0002 ±0.0031 (T(4S)), (14) 
-0.0068 ±0.0047 (D0). (15) 



Our combination of these values gives 



a", = ±0.0022 ± 0.0026. 



(16) 



The measurements of afj are performed by D0 pTj| and 
LHCb [111 collaborations: 



= -0.0104 ±0.0074 (DO), 



(17) 



-0.0024 ±0.0063 (LHCb). (18) 



Our combination of these values gives 
of, = -0.0058 ±0.0048. 



(19) 



Using the values (|16[> and (|19[) we obtain the allowed 
contributions to the inclusive muon and like-sign dimuon 

charge asymmetries from CPV in mixing of B° and B® 
mesons: 

os(fromo^) = (±0.008 ±0.009)%, (20) 

4s (from a£) = (±0.062 ±0.073)%, (21) 

a s (from< 1 ) = (-0.014 ± 0.012)%, (22) 

4s(fromafJ = (-0.111 ± 0.093)%, (23) 

and their sum 

a s (froma«[) = (-0.006 ±0.015)%, (24) 

4s (from a^) = (-0.049 ± 0.118)%. (25) 

For comparison, the SM prediction is 

4s (from aj in SM) = (-0.013 ± 0.002)%. (26) 
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In addition, the estimate of can be extracted from 
the measurement of CP violation in the B® — > J/ip(f> de- 
cay. The corresponding phase <f>1 cs is expected to change 
in the same way as the phase 4>\ 2 due to a new physics 
contribution^ @,[T3: 



Vl2 

<j>f s ! (SM) 



0f s (SM) + ^, 



(27) 

0f 2 (SM) + 0f 
■ sin(2^ s ) = -0.036 ± 0.002. (28) 

The current world average value of (jf s cs is Q : 



= -0.013 



+0.083 
-0.090- 



(29) 



From these expressions we get the following estimate of 
a|j from B° s — > J/ip(f> decay: 



aJ(from B° s -> J/ip<f>) = -0.0001 ± 0.0005. 



(30) 



However, we do not use the result (|3U|) because Eq. (f2"7| 
may be subject to large penguin corrections from new 
physics fl3j . 

Thus, the current experimental results constrain the 
contribution of CPV in mixing of B° and B® mesons to 
the measured like-sign dimuon charge asymmetry. 



B. CP violation in interference of B° decay with 
and without mixing 

The final states of the decay B°(B°) —> cedd are acces- 
sible from both B° and B°. Therefore, the interference of 
decays to these final states with and without B° mixing 
results in CPV Q. It turns out that this CPV produces 
a like-sign dimuon charge asymmetry. At the same time, 
its contribution is negligible in the inclusive muon charge 
asymmetry. 

To demonstrate this, let us consider an example of the 
process producing a positive dimuon pair 



PP 
B+ 
B° 



B+B°X, 
fj, + X, 

D+D-,D+ 



(31) 



The state /, = D + D~ is a CP-even eigenstate acces- 
sible from both B° and B° mesons. The D + fi + X 
decay produces a "wrong sign" muon, and contributes 
to the like-sign dimuon sample. The D~ yT X de- 
cay produces a "right sign" muon, and therefore does 
not contribute to the like-sign dimuon sample. The 
number of positive and negative muons from the decay 
B° —> D + D~ is the same, therefore there is no contribu- 
tion to the inclusive muon charge asymmetry from this 
decay. 

The decay rate of the meson that is initially produced 
as a B° is 



dr(S° -> ft 



dt 



tx exp (— T d t) x 



[1 + Si sin (Am d t) - d cos (Am d t)) . (32) 



The term proportional to Si is due to CPV in the in- 
terference of decays with and without mixing to D + D~ , 
of the meson that is initially produced as B°. The term 
proportional to Cj is due to the dirct CPV in B° decay. 
Neglecting the loop contributions to the decay amplitude, 
the coefficients Si and Ci in the SM are expressed as [l4| 



Si 
Ci 



-r\i sin(2/J), 



(33) 



0. 



where rji is the CP eigenvalue of the final state fi, and 
(3 is the angle of the Unitarity Triangle Q. The l oop 
diagrams can change this estimate by a few percent [15|, 
[3. For the B° D+D~ final state % = +1. If /, 
is not a CP eigenstate, the coefficient C, may be non- 
zero, but experimentally for the decays considered here 
Ci is negligible Q . Therefore, we omit it in the following 
discussion. Integrating (|32p we obtain the width of the 
decay to this final state 



rCB°->/i) cx 1 + 



SiXd 



..2 ' 



(34) 



T7' 



Now consider the CP-conjugate process producing a neg- 
ative dimuon pair 



PP 
B~ 
B° 



B~B°X, 



(,-X, 
D+D- 



H~X. 



(35) 



The decay rate of the meson that is initially produced as 
a B° is 



dT(B° -> f^ 



dt 



(X exp (-T d t) [1 - Si sin (Am d t)] , (36) 



and the partial width is 

T(B° -> fi) oc 1 



SjX d 

1 



X' 



(37) 



The like-sign dimuon charge asymmetry from this process 
is 

, x d 



A, 



1 + ^' 



(38) 



Numerically the absolute value of this asymmetry is 
large, because sin(2/3) = 0.679 ± 0.020 and x d = 0.770 ± 
0.008 @. 

Let us now obtain the contribution of the decay chan- 
nel fi with weight w^ d {fi) to the like-sign dimuon charge 
asymmetry As- This weight takes into account both the 
branching fraction of B° — s- fi —> fiX decay and detector- 
related efficiency of muon reconstruction. The weights 
for the various processes have been defined in Section [TTJ 
The probability that an initial b quark produces a "right 
sign" muon fi~ is [l[ 



P b cx 0.5w 3d {fi 



1 + Si- 



Xd 



- Wi a 

+W 2 b + 0.5(w 3s + W 3 f s +W i + W 5 ) 



(39) 
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The factor 0.5 in this expression corresponds to the state- 
ment that the number of positive and negative muons 
produced in the processes T3, T4 and T5 is the same. 
The probability that an initial b quark produces a "wrong 
sign" muon \i~ is 



P h oc 0.5w 3d (fi) 



Xd 



- w u 

+w 2a + 0.5(w 3s + w 3{s + W4 + w 5 ). 



(40) 



The number of events is proportional to PtPi- The 

number of events is obtained by replacing Si by 

— Si- The charge asymmetry from channels fi is then 



■r,i 



Pi 



l + xj 



PbPl 



(41) 



Thus, CPV in interference of B° decay with and with- 
out mixing produces a like-sign dimuon charge asymme- 
try, while it does not contribute to the inclusive muon 
charge asymmetry. The possible final states produced 
in the B°(B°) -s- cedd decay include D+D~, D*+D~, 
D+D*-, D*+D*~, J/t/'tt , J/tPn, J/ipp°, etc. For many 
of these final states the value Sj is measured experimen- 
tally [|[ and in all cases it is consistent with the SM 
value Si = — sin(2/3), which corresponds to the expecta- 
tion of the dominance of the CP-even final states in the 
B°(B°) -> cedd decay. 

On the contrary, the contribution of CPV in the 
b — > ccs transition producing the CP eigenstates, like 
the B° —> J/ipKs decay, should not contribute to the 
like-sign dimuon charge asymmetry, because for each 
CP-even final state there should be the correspond- 
ing CP-odd final state. For example, the contribution 
from the decay B° — > J/ipKs is canceled by the decay 
B° -> J/ij)K L . 

The weight w 3 d{fi) can be obtained by using the mea- 
sured branching fraction of B° — > fi decay. For example, 
the weight w 3 d{D + D~) can be computed using the fol- 
lowing expression 



W3d(fi) = fda 



Br(P° -> D+D-)2Bt(D+ -> fiX) 
Br(b -> fiX) ' 



(42) 



The coefficient w 3d {fi) by definition is normalised to the 
weight w\, which is proportional to Br(6 — > fJ-X). The 
value Br(6 — > /j,X) is the average branching fraction of 
the direct decay of all B hadrons to muon weighted with 
the relative production rate of different B hadrons at the 
hadron collider. Also by definition, the weight w 3 d{fi) 
includes both decays D + — > [i + X and D~ — > /i~X, hence 
the factors 0.5 in Eqs. ((39]) to ([Hjl. and the factor 2 in 
Eq. ((42]). The factor f d = 0.401 ±0.008 [§] is the fraction 
of B° plus B° in the admixture of 6-hadrons. The factor 
a is the ratio of detector acceptances of muons from D + 
and B° decay. Muons from D + and B° mesons have 
different detector acceptances because they have different 
kinematic distributions. 



Using the results of Ref. Q and other experimental 
values from Q we estimate the coefficient a from the 
following expressions: 



wi oc Br(6 
W3 oc Br(6 



ccX)Bv(ccqq' 



liX)a. (43) 



Here Br(6 — > ccX) is the branching fraction of B hadron 
decays producing cc pair. We use the experimental value 
Br(6 -> ccX) = 0.162 ± 0.032 which is obtained from 
Br(B ± AB°/S°/6-baryon mixture -> c/cX) = 1.162 ± 
0.032 Q. The quantity Bi(ccqq' — > fj,X) is the aver- 
age branching fraction of the direct decay of all charmed 
hadrons to muon weighted with the relative production 
rate of different pairs of c hadrons in B hadron decay. 
Using the values of corresponding branching fractions for 
the P ± /P°/Pg/&-baryon mixture from Ref. Q we ob- 
tain 



Br(6->/iX) = 0.107 ±0.003, 
Br(ccqq' -> fxX) = 0.164 ± 0.032. 

From these expressions we obtain 

Br(6 -> fiX) 



it's 



Br(6 -> ccX)Br{ccqq' -> pX) 
0.258 ±0.073. 



(44) 
(45) 



(46) 



This estimate of a does not take into account different 
kinematic distributions for the various decays c — > fiX. 
Therefore, a simulation of the D0 detector is required to 
obtain a more accurate value for a. 
Using Eqs. (g2} and (01]) we obtain 

Br(£?°->/j) Br(,/W^Q 
wM) = f d w 3 Qr{b - c _ x) Mc _ q _, - — } . (47) 

Using this observation, we estimate As by several 
methods. 

Estimate 1. Let us consider four measured decay 
channels 0: D+D~ with S = -0.87 ± 0.26, T./T = 
(2.11 ± 0.31) x 10" 4 ; P*(2010)+L>- with S = -0.61 ± 
0.19, r,/r = (6.1 ± 1.5) x 10~ 4 ; P*(2010)-P+ with S = 
-0.78±0.21, Ti/T w 6.1 x 10" 4 (our guess); and D*+D*~ 
with S = -0.76 ± 0.14, r,/r = (8.2 ± 0.9) x 10~ 4 . Using 
these numbers we obtain for the sum of these channels 

J2(SiBi{B° -> ftBitfi -> nX)) = 

i 

= -0.00044 ±0.00011. (48) 

Using this value and Eqs. (|4"T1) and (|47p. the contribution 
to As from these 4 channels is 



A s (4: channels) = (-0.028 ± 0.011)%. 



(49) 



Estimate 2. For this estimate we assume that the 
final states eded are mostly CP-even (rji = +1), which 
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is appropriate for £)(*)+!)(*) final states and confirmed 
by the experimental measurements. Then 



Br(P° -> cedd) 



(Br(6 ccs) + Br(6 -> ccd)) 



cd' 



(50) 



Eq. (|41j) becomes approximately 



A s 
S 



-0.5^/^ sin (2£)- 

Br(ccdJ -S- /zA") 
Br(ccgg' -> tzJT) ' 



■I',! 



PbPb 



'-6, 



(51) 



The factor 5 takes into account the fact that the final 
state cedd contains more mesons than the generic 
ccqq' state, and that the semileptonic branching fraction 
of Z? ± meson is about 2.7 times larger than that of all 
other charm hadrons. Using the known branching frac- 
tions of P-meson decays, we estimate S = 1.5±0.2. Using 
the values from Table U we obtain the following estimate 
of As from CPV in interference: 



A s = (-0.089 ±0.015)%. 



(52) 



This value gives an upper bound of the — As estimate, 
since in deriving it we assume that all cedd states are 
CP-even, which is clearly not the case. 

Estimate 3. In the SM the CP violation in mixing 
of neutral B mesons is small and the mass eigenstatcs 
of the B° system coincide with CP eigenstates. In that 
case [13 



Ar = AT CP = T(B° d ' even ) 



tv T}Q,odd\ 
1 i B d h 



(53) 



where T(B° d ' even ) (T(B° d ' odd )) is the width of P° decay 
to the CP-even (CP-odd) final states, respectively. As- 
suming that this difference is saturated by the P°(P°) — > 
cedd transition, we obtain the following estimate: 

£(Br(B° m) = 



sm(2f3)[Br(B° d even ) - Bi(B° d odd )} 
sin(2/3)Ar d /r d . 



(54) 



In the SM framework AT d /T d = (42 ± 8) x 10" 4 @. 
An alternative estimate of Ar^ uses the SM relation 
Ar d /Am d = Ar s /Am s given in Ref. @], page 1067. 
It results in a similar value ATd/Td = (44 ± 6) x 10~ 4 . 
Substituting the expression (JM)) in Eqs. (gTJ and ([i7]) 
we obtain 



A s = (-0.045 ±0.016)%. 



(55) 



Estimate 4. In this estimate we use the mea- 
sured value |Ar d /T d | = 0.015 ± 0.018 0. We replace 
E l (Br(P° -> fi)Si) by -sin(2/3)Ar d /r d and obtain 



A s = (-0.16 ±0.20)%, 



(56) 



All these estimates of As are consistent. In the fol- 
lowing, we use the estimate (|55|) . Comparing it with the 
experimental result ©, we conclude that CPV in inter- 
ference of B° decay with and without mixing accounts 
for a part of the observed like-sign dimuon charge asym- 
metry. The experimental uncertainty on Ar^ keeps open 
the possibility of a substantially larger contribution from 
this source. We also note that CPV in interference of B° 
is much larger than the SM prediction for CPV in mixing 
of B° and P° given in Eq. ||2BJ|. 

In Ref. \4 the like -sign dimuon charge asymmetry is 
measured in several sub-samples of events with an ad- 
ditional selection according to the muon impact param- 
eter. This selection effectively changes the contribution 
of muons coming from the oscillated B° decays. The es- 
timate (|55p after applying this selection is also modified 
and the contribution from CPV in interference can be en- 
hanced by selecting the dimuon events with large muon 
impact parameter. Therefore, the study of the depen- 
dence of the asymmetry As on the muon impact param- 
eter can provide an additional insight on this source of 
CPV. 



C. CPV in interference of B° decay with and 
without mixing 



Again we present several estimates. 

Estimate 1. Four channels of interest arc 
with iyr = (4.5 ± 1.4)%. [| The CPV parameters Si 
have not been measured, but in the SM should be ap- 
proximately 



Si 



sin(2/3 s ) 



-0.036. 



(57) 



Here, similarly to B° — > decays, we assume 

that these final states are mainly CP-even. For these 4 
decay channels we obtain 

A S (D^ + D^-) = (-0.0003 ± 0.0001(stat))%. (58) 

Estimate 2. We assume that the final states cscs 
are mostly CP-even (t^ = +1), which is appropriate for 
£)(*)+£>(*)- fi na i states . We take 



BrfP 



o 



cess) 



(Br(6 -> ccs) + Br(6 -> ccd)) 
Then Eq. (|41[) becomes approximately 



V. 



(59) 



-0.5w 3 f s V c 2 s sm{2f3 s )- 



1 + x* P b P h 
= (-0.0013 ± 0.0002(stat))%. 



(60) 



The absolute value of estimate (|60j) can be considered 
as an upper bound on the contribution to As from this 
source, because some of the P° —> cess final states are 
not CP-even. 
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Estimate 3. It is known that the four decay channels 
B s — > D s D s do not exhaust the contributions to 
Ar s Q. To obtain an estimate of the like-sign dimuon 
charge asymmetry from CPV in interference of B s we 
replace EjM^ -> fj)8 3 by - sin (2/3 s )Ar s /r s . We 
use the experimental value Ar s = 0.100 ± 0.013 ps _1 
and obtain 

A S {B° S ) = (-0.0009 ± 0.0003(stat))%, (61) 

which can be compared with @. 

In conclusion, CPV in interference of B° s decay with 
and without mixing is suppressed by the small values of 
sin (2f3 s ) and x s /(l + x 2 s ). 



D. Direct CPV in decay b 
followed by c — > 



acq (q = d or s), 
X. 



This type of CPV occurs due to the interference of 
the tree level and penguin diagrams with different strong 
phases and different weak phases. Let us consider, as 
an example, the decay B + — > D°D + , followed by D° — > 
[TX. Its branching fraction is (3.8 ± 0.4) x 10 -4 . The 
CP-violating asymmetry in this decay has not been mea- 
sured, but should be less than a — 0.1 because the pen- 
guin diagram is suppressed by one loop. The like-sign 
dimuon charge asymmetry from this channel, 



A s = af u aBr(B+ D°D+) 

Br(D° -> lTX)P b - Br(£H 



f i+X)P b 



P b P b Br(b -> fiX) 



(62) 



is less than 0.0002%. Considering also decays with q = s, 
which have a CP- violating asymmetry suppressed by one 
loop x A 2 , we conclude that the direct CP violation in the 
decays measured so far have a negligible contribution to 
the like-sign dimuon charge asymmetry. 



TABLE II: Contributions to As allowed by experiments. 





Process 


allowed As 


A 


Mixing of B" 


(+0.062 ± 0.073)% 


A 


Mixing of B° 


(-0.111 ±0.093)% 


B 


Interference of B° 


(-0.045 ±0.016)% 


C 


Interference of B° 


(-0.0009 ± 0.0003)% 


D 


CPV in b — >• ccq decays 


(±0.000 ± 0.001)% 


E 


0(j,) in b — > fiX decays 


(-0.17 ±0.43)% 


E 


<j( c ) in c — > ^X decays 


(-0.07 ±0.19)% 



contribution of and a( c ) to as and As are 

toiO(i,) ± (w 2 + w 3 +w 5 + io 6 )a( c ) 



wi ± w 2 ± w 3 ± Wi ± 1V5 ± we 
= 0.535a (b) ± 0.454a (c ), 



(63) 



and 



WiaP b ± WibP b 



(6) 



w 2b P- b ± w 2a P b ± 0.5(w 3 ± w 5 )(P- b + P b ) 



P b P b 



1.424a 



(6) 



0.525a 



(c) 



-a(c) 
(64) 



If the only asymmetry is ar b \, then As /as — 2.66. If 
the only asymmetry is ai c \, then As /as = 1.16. Taking 
as from (p}, we obtain the following estimates for the 
contribution of direct CPV to As: 

,4 S (from a (b) ) = (-0.17 ± 0.43)%, 

,4 S (from a (c) ) = (-0.07 ± 0.19)%. (65) 

Thus, the "closure test" (HJ begins to constrain the 
contributions of a( c ) to the like-sign dimuon charge asym- 
metry As- 



IV. CONCLUSIONS 



E. Direct CPV in semileptonic decays of b and c 
quarks. 

In this subsection we assume that the like-sign dimuon 
charge asymmetry is due solely to CPV in direct semilep- 
tonic decays of charged and neutral hadrons containing 
b or c quarks. This type of CP violation vanishes in the 
lowest order due to the CPT symmetry [l^. The sec- 
ond order calculations give extremely small value for the 
asymmetry As of the order of 10~ 9 [l8|. Despite such a 
strong theoretical constraint, the possibility of large con- 
tribution from this source is discussed in Ref . [l9[ . There 
is no direct experimental measurements of this CPV in 
semileptonic B decays, and the only experimental limi- 
tation can be derived from Eqs. ([T]) and @. Let a/ b \ 
(ar c \) be the flavor averaged CP violating charge asym- 
metry of direct semileptonic decays of b (c) quarks. The 



We have considered several possible causes of the mea- 
sured like-sign dimuon charge asymmetry As, and ob- 
tained their experimental constraints. A summary is 
presented in Table [Til We find that standard model CP 
violation in interference of decays with and without mix- 
ing of B° to flavor non-specific states /,-, followed by the 
decay ft —> \xX, contributes 

A S (B°) = (-0.045 ± 0.016(stat))% (66) 

to the like-sign dimuon charge asymmetry As- CP vio- 
lation in interference does not contribute to the inclusive 
muon charge asymmetry and therefore is compatible with 
the observation that as is consistent with zero. 

Among all other possible sources of the dimuon charge 
asymmetry, only the direct CP violation in semileptonic 
b- and c-hadron decays is not yet limited experimentally. 
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It is very small in the SM, but, until experimentally mea- 
sured, this source of dimuon charge asymmetry cannot be 
excluded. Our estimate of this source is derived from the 
D0 measurements ([T]) and The exclusive measure- 
ment of this type of CP violation is required to improve 
this constraint. 

Taking into account the additional SM source of 
dimuon charge asymmetry (|66[) identified in this paper, 
the combination of D0 measurements ([T]) and ([2j be- 
comes consistent with the SM expectation within 3 stan- 
dard deviations. Still the difference between ©, and 
(Hp and HBBJI, (-0.32 ± 0.14 (tot))%, leaves some room 



for new physics CP violation in B° and B® mixing, in 
the interference of B° and B® decays with and without 
mixing, or in scmilcptonic decays of b and c hadrons. A 
deviation in the value of AI^ from the SM prediction 
could also contribute to the difference between the ob- 
served and expected like-sign dimuon charge asymmetry. 

As DO collaborators deeply involved in the dimuon 
asymmetry analysis, we thank the D0 Collaboration for 
inspiring, commenting, and supporting this work. We 
thank Michael Gronau and Jonathan L. Rosner for very 
useful comments and suggestions to our paper. 
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